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Abstract 

Cross sections for e~p charged current deep inelastic scattering have been mea- 
sured at a centre-of-mass energy of 318 GeV with an integrated luminosity 
of 16.4 pb" 1 using the ZEUS detector at HERA. Differential cross-sections 
da/dQ 2 , da/dx and da/dy are presented for Q 2 > 200 GeV 2 . In addition, 
d 2 a/dxdQ 2 was measured in the kinematic range 280 GeV 2 < Q 2 < 30 000 GeV 2 
and 0.015 < x < 0.42. The predictions of the Standard Model agree well with 
the measured cross sections. The mass of the W boson, determined from a fit to 
da/dQ 2 , is M w = 80.3 ± 2.1 (stat.) ± 1.2 (syst.) ± 1.0 (PDF) GeV. 
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1 Introduction 



Deep inelastic scattering (DIS) of leptons on nucleons has been vital in the development 
of our understanding of the structure of the nucleon. In the Standard Model (SM), 
charged current (CC) DIS is mediated by the exchange of the W boson. In contrast 
to neutral current (NC) interactions, where all quark and antiquark flavours contribute, 
only up-type quarks and down-type antiquarks participate at leading order in e~p CC DIS 
reactions. This makes such interactions a powerful tool for flavour-specific investigation 
of the parton distribution functions (PDFs). Since only left-handed quarks and right- 
handed antiquarks contribute to CC DIS at HERA, the distribution of the electron-quark 
centre-of-mass scattering angle, 8*, is a sensitive probe of the chiral structure of the weak 
interaction. 

Measurements of the CC DIS cross sections at HERA have been reported previously by 
the HI [1,2] and ZEUS [3,4] collaborations. These data extended the kinematic region 
covered by fixed-target neutrino-nucleus scattering experiments [5] by about two orders 
of magnitude in the negative square of the four-momentum transfer, Q 2 . In addition, the 
double-differential e + p CC DIS cross section, d 2 a/dxdQ 2 , where x is the Bjorken scaling 
variable, was measured for the first time at high Q 2 by the HERA collider experiments 
[6,7]. The mass of the exchanged boson in the space-like domain, extracted from a fit 
to the differential cross-section da/dQ 2 , was consistent with the mass of the W boson 
measured in time-like processes at LEP and at the Tevatron [8]. 

This paper presents measurements of the e~p CC DIS single-differential cross-sections 
da/dQ 2 , da/dx and da/dy, as well as d 2 a/dxdQ 2 . The results are compared to the 
expectations of the SM. The measurements are based on 16.4 pb _1 of data collected 
during the running periods in 1998 and 1999 when HERA collided electrons of energy 
27.5 GeV with protons of energy 920 GeV, yielding a centre-of-mass energy of 318 GeV. 
The data represent an increase of a factor of 20 in integrated luminosity over the previous 
ZEUS e~p measurement [4]. Cross sections for e~p CC DIS were reported recently by the 
HI collaboration [9]. 



2 Standard Model prediction 

The electroweak Born-level CC DIS differential cross section, d 2 a^ rn /dxdQ 2 , for the re- 
action e~p — > v e X, with longitudinally unpolarised beams, can be expressed as [10] 



d 2 cr^ rn (e-p) _ G 2 F 

dxdQ 2 " Awx (Q 2 + M^ 2 



[Y + F 2 cc (x, Q 2 ) - y 2 F? c (x, Q 2 ) + Y.xF 3 cc (x, Q 2 )] (1) 



1 



where Gp is the Fermi constant, M\y is the mass of the W boson, x is the Bjorken scaling 
variable, y = Q 2 /xs and Y± — 1 ± (1 — y) 2 . The centre-of-mass energy in the electron- 
proton collision is given by s/s = 2y/E e E p , where E e and E p are the electron and proton 
beam energies, respectively. The inelasticity y is related to 9* by y = 1/2(1 — cos#*). The 
structure functions F^ and xF^ , at leading order in QCD, may be written in terms 
of sums and differences of quark and antiquark PDFs. For longitudinally unpolarised 
beams, 

F 2 CC = x[u(x, Q 2 ) + c(x, Q 2 ) + d(x, Q 2 ) + s(x, Q% 

xF^ c = x[u(x, Q 2 ) + c{x, Q 2 ) - d(x, Q 2 ) - s(x, Q% 

where, for example, the PDF u(x, Q 2 ) gives the number density of up quarks with mo- 
mentum fraction x at a given Q 2 . The longitudinal structure function, F^ c , is zero at 
leading order in QCD. At next-to-leading-order (NLO) in QCD F^ c is non-zero but gives 
a negligible contribution to the cross section except at values of y close to 1, where it can 
be as large as 10%. Since the top-quark mass is large and the off-diagonal elements of 
the CKM matrix are small [8], the contribution from the third-generation quarks may be 
ignored in CC DIS at HERA [11]. Since the w-quark density, which is well constrained 
by NC DIS data, dominates in e~p CC DIS, the uncertainties coming from the PDFs are 
small. The uncertainty in the prediction for da/dQ 2 , for example, increases from 2% to 
5% over the Q 2 range of the present measurements. In the following, uncertainties in the 
predicted cross sections were obtained using the ZEUS NLO QCD fit [12]. 

3 The ZEUS experiment 

A detailed description of the ZEUS detector can be found elsewhere [13]. A brief outline 
of the components most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [14], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organised in nine superlayers covering the polar- 
angleQ region 15° < 9 < 164°. The relative transverse-momentum resolution for full-length 
tracks is a(p T ) jp T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. The position of the 
interaction vertex along the beam direction can be reconstructed from the CTD tracks 
with a resolution of about 1 cm in CC events. 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The polar angle, 9, is measured with respect to the proton beam direction. The 
coordinate origin is at the nominal interaction point. 
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The high-resolution uranium-scintillator calorimeter (CAL) [15] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL relative 
energy resolutions, as measured under test-beam conditions, are a(E)/E = 0.18/ VE for 
electrons and a(E)/E = 0.35/y/E for hadrons (E in GeV). The timing resolution of the 
CAL is better than 1 ns for energy deposits greater than 4.5 GeV. The position of the 
interaction vertex along the beam direction can be reconstructed from the arrival time 
of energy deposits in FCAL. The resolution is about 9 cm for events with FCAL energy 
above 25 GeV, improving to about 7 cm for FCAL energy above 100 GeV. 

An instrumented-iron backing calorimeter (BAC) [16] surrounds the CAL and can be used 
to measure energy leakage and identify muons. Muon chambers in the forward, barrel and 
rear [17] regions are used in this analysis to veto background events induced by cosmic-ray 
or beam-halo muons. 

The luminosity was measured using the Bethe-Heitler reaction ep — > e^p. The photons 
were measured by the luminosity monitor [18], a lead-scintillator calorimeter placed in 
the HERA tunnel 107 m from the interaction point in the electron beam direction. 

4 Monte Carlo simulation 

Monte Carlo (MC) simulation was used to determine the efficiency for selecting events 
and the accuracy of kinematic reconstruction, to estimate the ep background rates and 
to extract cross sections for the full kinematic region. A sufficient number of events was 
generated to ensure that the statistical uncertainties arising from the MC simulation were 
negligible compared to those of the data. The MC samples were normalised to the total 
integrated luminosity of the data. 

The ZEUS detector response was simulated with a program based on geant 3.13 [19]. 
The generated events were passed through the simulated detector, subjected to the same 
trigger requirements as the data, and processed by the same reconstruction programs. 

Charged current DIS events, including electroweak radiative effects, were simulated using 
the HERACLES 4.6.1 [20] program with the djangoh 1.1 [21] interface to the MC gener- 
ators that provide the hadronisation. Corrections for initial-state radiation, vertex and 
propagator corrections and two-boson exchange are included in HERACLES. The mass of 
the W boson was calculated using the PDG [8] values for the fine structure constant, 
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the Fermi constant, the mass of the Z boson and the mass of the top quark, and with 
the Higgs-boson mass set to 100 GeV. The colour-dipole model of ARIADNE 4.10 [22] 
was used to simulate the 0{ots) plus leading logarithmic corrections to the result of the 
quark-parton model. As a systematic check, the meps model of lepto 6.5 [23] was used. 
Both programs use the Lund string model of jetset 7.4 [24] for the hadronisation. A 
set of NC events generated with djangoh was used to estimate the NC contamination 
in the CC sample. Photoproduction background was estimated using events simulated 
with herwig 5.9 [25]. The background from W production was estimated using the 
epvec [26] generator, and the background from production of charged-lepton pairs was 
generated with the lpair [27] program. 

5 Reconstruction of kinematic variables 

The principal signature of CC DIS at HERA is the presence of a large missing transverse 
momentum, Pr.miss, arising from the energetic final-state neutrino that escapes detection. 
The quantity -Fr.miss was calculated from 



where the sums run over all calorimeter energy deposits, Ei (uncorrected in the trigger, 
but corrected in the offline analysis for energy loss in inactive material etc. [28]) and 0j 
and <f>i are the polar and azimuthal angles of the calorimeter deposits as viewed from the 
interaction vertex. The polar angle of the hadronic system, 7^, is defined by cos 7^ = 



(^r,miss - $ 2 )/(PtMs S + where 6 = E(Ei ~ EiCOsOi) = - P z )i. In the naive 



quark-parton model, 7^ is the angle through which the struck quark is scattered. Finally, 
the total transverse energy, Et, is given by Et = ^2,Ei sm ®%- 

The kinematic variables were reconstructed using the Jacquet-Blondel method [29]. The 
estimators of y, Q 2 and x are: y JB = 5/{2E e ), Q 2 JB = ir )miss /(l - Vjb), and x JB = 
Q 2 j B /(syjB)- 

6 Event selection 

Charged current DIS candidates were selected by requiring a large -Pr.miss and a recon- 
structed event vertex consistent with an ep interaction. The main sources of background 
come from NC scattering and high-i?T photoproduction. The energy resolution of the 




Ei sin 6i cos (f>i I + / Ei sin 6{ sin 0; 
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CAL or the energy that escapes detection can lead to significant missing transverse mo- 
mentum. Events not from ep collisions, such as beam-gas interactions, beam-halo muons 
or cosmic rays can also cause substantial apparent imbalance in the transverse momentum 
and so constitute other sources of background. The selection criteria described below were 
imposed to separate CC events from all backgrounds. 

When the current jet lies in the central region of the detector, i.e. jh is large, tracks in the 
CTD can be used to reconstruct an event vertex, strongly suppressing non-ep backgrounds. 



The procedure designed to select these events is described in Section |6.2| . For CC events 
with small 7^, the charged particles from the hadronic final state are often outside the 
acceptance of the CTD. Such events populate the high-x region of the kinematic plane. 
The algorithm designed specifically to select such events is described in Section |0T3| . The 
events were classified first according to 70, the value of 7^ measured with respect to the 
nominal interaction point. Subsequently, the kinematic quantities were recalculated using 
the Z-coordinate of the event vertex (Zvtx) determined from either CTD tracks or the 
calorimeter-timing information discussed in Section |3|. 

In the data-taking period, the HERA proton beam contained a significant number of off- 
axis protons. Such particles passed through the ZEUS interaction region displaced hori- 
zontally by several millimeters from the nominal trajectory and gave rise to an especially 
pernicious class of off-axis beam-gas interactions. In order to remove this background, 
several selection requirements had to be made more restrictive than was the case in the 
previous ZEUS study [6]. 



6.1 Trigger selection 

ZEUS has a three-level trigger system [13]. At the first level, events were selected using 
criteria based on the energy, transverse energy and missing transverse momentum mea- 
sured in the calorimeter. Generally, events were triggered with a low threshold on these 
quantities when a coincidence with CTD tracks from the event vertex was required, while 
a higher threshold was necessary for events with no CTD tracks. The latter class of events 
has a hadronic final state boosted forward outside the CTD acceptance. Typical threshold 
values were 5 GeV (8 GeV) in missing transverse momentum, or 11.5 GeV (21 GeV) in 
transverse energy for events with (without) CTD tracks. 

At the second level, timing information from the calorimeter was used to reject events 
inconsistent with the bunch-crossing time. In addition, the topology of the CAL energy 
deposits was used to reject background events. In particular, since the resolution of the 
missing transverse momentum is better at the second level than at the first level, a tighter 
cut of 6 GeV (9 GeV for events without CTD tracks) was made. 
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At the third level, track reconstruction and vertex finding were performed and used to 
reject candidate events with a vertex inconsistent with the beam interaction envelope. 
Cuts were applied to calorimeter quantities and reconstructed tracks to reduce beam-gas 
contamination further. 

6.2 Offline selection based on a CTD vertex 

Events with 7 > 0.4 rad were required to have a vertex reconstructed from CTD tracks 
and to satisfy all of the following criteria: 

• |Zvtx| < 50 cm 

the primary vertex reconstructed from the CTD tracks was required to be within the 
range consistent with the ep interaction region; 

• P T , miss > 12 GeV and miss > 10 GeV 

-^Tmiss i s the missing transverse momentum calculated excluding the FCAL towers 
closest to the beam hole. The Pf miss cu ^ strongly suppresses beam-gas events while 
maintaining high efficiency for CC events; 

• Tracking requirement 

at least one track associated with the event vertex must have transverse momentum 
in excess of 0.2 GeV and a polar angle in the range 15° to 164°. In order to remove 
off-axis beam-gas background, a cut was also applied in two dimensions on the number 
of such "good" tracks, N goo< i, versus the total number of tracks, [30]; 

• Rejection of photoproduction 

Pt, miss / ' Et > 0.5 was required for events with Pt, miss < 30 GeV. This cut selected a 
collimated energy flow, as expected from a single scattered quark. No Pt,-hus&/Et re- 
quirement was imposed on the events with Pr.miss > 30 GeV. In addition, the difference 
between the direction of the (Px,Py) vector calculated using CTD tracks and that 
obtained using the CAL was required to be less than 0.5 radians for Pr, m iss < 30 GeV 
and less than 2.0 radians for Pr, m iss > 30 GeV; 

• Rejection of NC DIS 

NC DIS events in which the scattered electron or jet energy was poorly measured 
could have a large apparent missing transverse momentum. To identify such events, 
a search for candidate electrons was made using an electron-finding algorithm that 
selects isolated electromagnetic clusters in the CAL [31]. Candidate electron clus- 
ters within the CTD acceptance were required to have an energy above 4 GeV and a 
matching track with momentum larger than 25% of the cluster energy. Clusters with 
6 > 164° were required to have a transverse momentum exceeding 2 GeV. Events with 
a candidate electron satisfying the above criteria and 5 > 20 GeV were rejected. For 
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contained NC events, 5 peaks at 2E e = 55 GeV. This cut was only applied for events 
with Pr.miss < 30 GeV; 

• Rejection of non-ep background 

beam-gas events typically gave CAL times that were inconsistent with the bunch- 
crossing time. Such events were rejected. Muon-finding algorithms based on CAL 
energy deposits or muon-chamber signals were used to reject events produced by cosmic 
rays or muons in the beam halo. 



6.3 Offline selection based on a CAL vertex 

For events with 70 < 0.4 rad, the following criteria were imposed: 

• |Zvtx| < 50 cm 

Zvtx was reconstructed from the measured arrival time of energy deposits in FCAL [32] ; 

• Ptmss > 25 GeV and P^ miss > 25 GeV 

to remove off-axis beam-gas events, the conditions on missing transverse momentum 
were tightened, compensating for the relaxation of the requirements on tracks; 

• Rejection of non-ep background 



the timing and muon-rejection cuts described in Section |6.2| were used. A class of 
background events that were especially troublesome in this selection branch arose from 
beam-halo muons that produced a shower inside the FCAL. To reduce this background, 
topological cuts on the transverse and longitudinal shower shape were imposed; these 
cuts rejected events in which the energy deposits were much more strongly collimated 
than for typical hadronic jets. 



6.4 Final event sample 

To restrict the sample to regions where the resolution in the kinematic variables was 
good and the backgrounds small, the requirements Q 2 JB > 200 GeV 2 and yjs < 0.9 
were imposed. All events were visually inspected, and nine cosmic-ray and halo-muon 
events were removed. The final e~p sample consisted of 650 events, to be compared with 
660 predicted by the MC simulation. 

The resolution in Q 2 is ~ 25% over the entire Q 2 range. The resolution in x improves 
from ~ 25% at x = 0.01 to ~ 10% at x = 0.5. The resolution in y is ~ 13% over the 
entire range. 

Figure [l] compares the distributions of data events entering the final CC sample with 
the MC expectation for the sum of the CC signal and ep background events. The MC 
simulation gives a good description of the data. 
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7 Cross-section determination and systematic uncer- 
tainties 



7.1 Cross-section determination 

Monte Carlo events were generated according to Eq. flU), including electroweak radiative 
effects. The value of the cross section, at a fixed point within a bin, was obtained from 
the ratio of the number of observed events, from which the estimated background had 
been subtracted, to the number of events predicted by the MC simulation, multiplied by 
the cross section obtained using Eq. (|l|). Consequently, the acceptance, bin-centring and 
radiative corrections were all taken from the MC simulation. 

7.2 Systematic uncertainties 

The major sources of systematic uncertainty in the cross sections are discussed below: 

• Uncertainty of the calorimeter energy scale 

the hadronic energy scale and the associated uncertainty were determined, using NC 
DIS events, from the ratios of the total hadronic transverse momentum, Pr,h&d, to 
Pt,da and Pr, e , where Pt,ba = yj Qba 0- ~~ 2/da) is the transverse momentum obtained 
from the double-angle method [33] and P-r,e is the measured transverse momentum of 
the scattered electron. In order to restrict the hadronic activity to particular polar- 
angle regions, a sample of NC DIS events with a single hadronic jet was selected. By 
applying suitable cuts on the location of the current jet and evaluating Pt, had/ Pt,t> a 
and Pr,had/ PT,e event by event, the hadronic energy scales of the FCAL and BCAL 
were determined. The responses of the HAC and EMC sections of the individual 
calorimeters were determined by plotting Pt, h&d / Pt, da and Pt, had / Pt,c as a function 
of the fraction of the hadronic energy measured in the EMC section of the calorimeter. 
In each case, the uncertainty was found by comparing the determinations from data 
and MC. In order to study the hadronic energy scale in the RCAL, a sample of 
diffractive DIS events was selected. Such events are characterised by a large gap in the 
hadronic energy flow between the proton remnant and the current jet. PT,had/ Pt,da 
was evaluated event-by-event for events with hadronic activity exclusively in the RCAL 
and the energy scale and associated uncertainty determined. 

The relative uncertainty in the hadronic energy scale was determined to be 2% for the 
RCAL and 1% for the FCAL and BCAL. Varying the energy scale of the calorime- 
ter sections by these amounts in the detector simulation induces small shifts of the 
Jacquet-Blondel estimators for the kinematic variables. Varying the energy scale of 



S 



each of the calorimeters simultaneously up or down by these amounts in the detector 
simulation gives one estimate of the systematic uncertainty A second estimate was 
obtained by increasing (decreasing) the FCAL and RCAL energy scales together while 
the BCAL energy scale was decreased (increased). A third estimate was made by 
simultaneously increasing the energy measured in the EMC section of the calorimeter 
by 2% and decreasing the energy measured in the HAC section by 2% and vice-versa. 
This was done separately for each of the calorimeters. The final systematic error at- 
tributed to the uncertainty in the hadronic energy scale was obtained by taking the 
quadratic sum of these three estimates. The resulting systematic shifts in the mea- 
sured cross sections were typically within ±5%, but increase to +10% in the highest 
y bin and ±(10 — 20)% in the highest Q 2 and x bins; 

• Energy leakage 

four percent of the accepted events have a measurable energy leakage from the CAL 
into the BAC. The average transverse-energy leakage for these events is 5% of the 
Et observed in the CAL. Both the fraction of events with leakage and the average 
leakage are well modelled by the MC simulation and the effect on the cross-section 
measurement is negligible; 

• Variation of selection thresholds 

for the majority of the kinematic bins, varying the selection thresholds by ±10% in 
both data and MC resulted in small changes in the measured cross sections. However, 
varying the Pp, miss / E T threshold gave a change in the cross section of around ±10% 
at low Q 2 in da/dQ 2 and up to ±(10 — 20)% in the bins of the double-differential 
cross section. Varying the Pf m i SS cu ^ S ave a change of approximately ±10% in the 
single-differential cross-section da/dy at the lowest y. Varying the Pr,miss cut gave 
a change of —12% in the lowest-Q 2 bin of da/dQ 2 . The tracking requirement for 
Pr,miss < 30 GeV was tightened and gave changes in the cross section of ±8% and 
— 12% in the lowest-Q 2 bins of da/dQ 2 and —10% in the highest y bin of da/dy. 
Changes of ±(10 — 20)% were observed in the lowest-Q 2 bins of the double-differential 
cross section; 

• Uncertainty in the parton-shower scheme 

the meps model of lepto was used to calculate the acceptance instead of the ARIADNE 
model. The largest effects were observed in the bins of low Q 2 (±8%), highest and 
lowest y (±5%) and low x (±4%). The largest effect in the double-differential cross 
section was seen in the low-Q 2 and low-x bins, where it amounted to ±(10 — 13)%; 

• Background subtraction 

the uncertainty in the photoproduction background was estimated by fitting a linear 
combination of the Pr, ™;^ / Et distributions of the signal and the background MC 
samples to the corresponding distribution in the data, allowing the normalisation 
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(Nphp) of the photoproduction MC events to vary. No cut on P T> m \ ^/E T was applied 
for this check. Nphp was varied in the range Np/jp/2 to 2-N php. This range corresponds 
to the uncertainty given by the fit and changes in the cross sections were found to be 
within ±1%; 

• Trigger efficiency 

the simulation of the efficiency of the first-level trigger as a function of -P^miss was 
examined by using data events triggered by an independent trigger branch that was 
highly efficient for CC events, since it was based on CAL energy sums. The difference 
between the efficiencies calculated from the data and from MC events had a negligible 
effect on the measured cross section; 

• Choice of parton distribution functions 

the CC MC events were generated using the CTEQ5D PDFs [34]. The ZEUS NLO 
QCD fit [12] was used to examine the influence of variations of the PDFs on the cross- 
section measurement. Monte Carlo events were re- weighted to reflect the uncertainty 
from the fit and new acceptance-correction factors were computed. The change in the 
measured cross section was typically < 1%, except at low Q 2 where it was —2% and 
at high x where it was +3%; 

• The effect of F L 

the djangoh program neglects the Fl contribution to d 2 a/dxdQ 2 when generating CC 
events. The corresponding effect on the acceptance-correction factors was evaluated 
by re- weighting MC events with the ratio of the cross sections with and without F^. 
The largest effect was observed in the highest-?/ bin where it was —2%; 

• Uncertainty in the radiative correction 

the magnitude of the 0{a) electroweak corrections to CC DIS have been discussed by 
several authors [35,36]. Various theoretical approximations and computer codes gave 
differences in the CC cross sections of typically < ± (1 — 2)%. The differences can be 
as large as ±(3 — 8)% at high x and at high y. This uncertainty was not included in 
the total systematic uncertainty. 



The individual uncertainties were added in quadrature separately for the positive and 
negative deviations from the nominal cross-section values to obtain the total systematic 
uncertainty listed in Tables [l] and |[ The uncertainty on the luminosity of 1.8% was not 
included. 
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8 Results 



The single-differential cross-sections da/dQ 2 , da/dx and da/dy for Q 2 > 200 GeV 2 are 
shown in Fig. |2| and compiled^] in Table [|. The cross sections da/dQ 2 and da/dx were 
extrapolated to the full y range using the SM predictions with CTEQ5D PDFs. The SM 
cross sections derived from Eq. (0) using the ZEUS NLO QCD fit, CTEQ5D [34] and 
MRST(99) [37] parameterisations of the PDFs are also shown, together with the ratios of 
the measured cross sections to the SM cross section evaluated with the ZEUS NLO QCD 
fit. The Standard Model gives a good description of the data. 

The reduced double-differential cross section, a, is defined by 

d 2 a 
dx dQ 2 

At leading order in QCD, a(e~p — > v e X) depends on the quark momentum distributions 
as follows: 

a(e~p -> v e X) = x[u + c+ (I- y) 2 (d + s)] . (2) 

The reduced cross sections are displayed versus Q 2 and x in Figs. ||] and respectively, 
and compiled 2 in Table @. The predictions of Eq. ([!]), evaluated using the ZEUS NLO 
QCD fit, give a good description of the data. The PDF combinations (u + c) and (d+ s), 
obtained in the MS scheme from the ZEUS NLO QCD fit, are shown separately in Fig. 

The W boson couples only to left-handed fermions and right-handed antifermions. Thus, 
in the quark-parton- model description of e~p CC DIS, the Z component of the total 
angular momentum of the initial-state electron-quark system is zero for e~q scattering 
and one unit of angular momentum for e~q scattering. Therefore, the angular distribution 
of the scattered quark in e~q CC DIS will be flat, while it will exhibit a (1 + cos#*) 2 
distribution in e~q scattering. Since (1 — y) 2 = 1/4(1 + cos^*) 2 the helicity structure of 
CC interactions can be illustrated by plotting the reduced double-differential cross section 
of Eq. (Q) versus (1 — y) 2 in bins of fixed x. At leading order in QCD, in the region of 
approximate scaling, this yields a straight line, the intercept of which gives the (u + c) 
contribution, while the slope gives the (d + s) contribution. 

Figure |5] shows the e~p CC DIS data, compared to the previously published e + p data [6]. 
At large x, e~p CC DIS is sensitive to the w-valence-quark PDF while e + p CC DIS is 
sensitive to the ci-valence-quark PDF. At leading order in QCD, the reduced double- 
differential cross section for e + p CC DIS can be written as 

a(e + p -> v e X) =x[u + c+(l- y) 2 (d + s)] , 

2 Tables || and || contain details of the systematic uncertainties that are correlated between cross-section 
bins. 



G 2 
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permitting a similar interpretation for the intercept and slope in terms of the appropriate 
parton densities. Scaling violations can be observed in the theoretical prediction as (1— y) 2 
approaches 1. The data agree with the expectation of the SM from the ZEUS NLO QCD 
fit, which is also shown. 

The total cross section for e~p CC DIS in the kinematic region Q 2 > 200 GeV 2 is 
a^ T (Q 2 > 200 GeV 2 ) = 66.7 ± 2.7(stat.)lJ;?(syst.) pb. 

The result is in good agreement with the SM expectation, evaluated using the ZEUS NLO 
QCD fit, of 69.011:1 P b - 

9 Electroweak analysis 

Equation ([!]) shows that the magnitude of the CC DIS cross section is determined by 
Gp and the PDFs. The fall in the cross section with increasing Q 2 is dominated by the 
propagator term, M^r/(Q 2 + M^) 2 . Figure |] shows that the Q 2 dependence of the PDFs 
is small by comparison. An electroweak analysis, performed by fitting da/dQ 2 with Gp 
fixed at the PDG [8] value of 1.16639 ■ 10~ 5 GeV" 2 and M w treated free parameter 
gives 

M w = 80.3 ± 2.1 (stat.) ± 1.2 (syst.) ± 1.0 (PDF) GeV, 

where the third uncertainty was estimated by varying the PDFs within the uncertain- 
ties given by the ZEUS NLO QCD fit [12]. The systematic uncertainty includes con- 
tributions from the sources identified in Section [T2] and the uncertainty on the mea- 
sured luminosity. This measurement, in the space-like region, is in good agreement 
with the more precise measurements of jy-boson production in the time-like region [8]. 
This measurement is also in good agreement with the previous ZEUS measurement of 
M w = 81.4l 2 ;£(stat.) ± 2.0(syst.)l3;o(PDF) GeV obtained by a fit to da/dQ 2 for e + p CC 
DIS at y/s = 300 GeV [6] . The present measurement has a smaller statistical uncertainty 
than the previous result, despite a lower integrated luminosity, due to the larger cross 
section for CC DIS in e~p scattering compared to e + p scattering. The measurement also 
benefits from improved understanding of the detector, leading to a smaller systematic 
uncertainty. Since the w-quark density, which is well constrained by NC DIS data, domi- 
nates in e~p CC DIS, the uncertainty coming from the PDFs is smaller than in the case 
of e + p CC DIS, where the <i-quark density dominates. 
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10 Summary 



Differential cross sections for charged current deep inelastic scattering, e~p — > iy e X, have 
been measured for Q 2 > 200 GeV 2 using 16.4 pb _1 of data collected with the ZEUS detec- 
tor during the period 1998 to 1999. The double-differential cross-section d 2 a/dx dQ 2 has 
been measured in the kinematic range 280 GeV 2 < Q 2 < 30 000 GeV 2 and 0.015 < x < 0.42. 
The chiral structure of the Standard Model was investigated by plotting the double- 
differential cross section as a function of (1— y) 2 . The Standard Model gives a good descrip- 
tion of the data. The mass of the W boson has been determined, from the Q 2 dependence 
of the measured cross section, to be M w = 80.3 ± 2.1 (stat.) ± 1.2 (syst.) ± 1.0 (PDF) GeV 
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Table 1: Values of the differential cross-sections da/dQ 2 , da/dx and da/dy. 
The following quantities are given: the range of the measurement; the value at 
which the cross section is quoted; the number of data events, -/Vdata/ the number 
of expected background events, Nbg and the measured cross section, with statistical 
and systematic uncertainties. 



17 



Q 2 ( GeV 2 ) 


X 


Ndata 


N BG 




280 





032 


7 





3 


n fu+ - 36 +°- 13 

u - u ^-0.25 -0.06 


530 





015 


10 








i i q+0.49 +0.14 
i - io -0.36 -0.26 


530 





032 


12 





1 


n qi +0.35 +0.08 
u - Ji -0.26 -0.16 


530 





068 


8 








n oq+0.41 +0.14 
u.oo_ 2g _ 19 


950 





015 


17 





7 


1 


51 ±0 


oq+0.25 
oy -0.28 


950 





032 


16 





1 





75 ±0 


1 q+0-08 

iy -0.08 


950 





068 


29 





2 


1 


01 ±0 


1O+0.13 
iy -0.13 


950 





130 


17 





1 





79 ±0 


10+0.11 

iy -0.09 


1700 





032 


31 





2 





87 ±0 


I6+0.11 

iu -0.05 


1700 





068 


51 





3 





81 ±0 


1 9+O.O2 
iz -0.05 


1700 





130 


32 





1 





75 ±0 


1 q+0.03 

io -0.05 


1700 





240 


18 











47 ±0 


1 i+O.Ol 
1 -1-0.02 


3000 





068 


56 





5 





71 ±0 


1 n+0.04 
iu -0.01 


3000 





130 


34 





1 





58 ±0 


1 n+0.03 

iu -0.02 


3000 





240 


27 











49 ±0 


QO+0.01 
ua -0.01 


3000 





420 


8 








n 91 +0.10 +0.02 
u - zi -0.07 -0.02 


5300 





068 


31 





3 





53 ±0 


1 n+o.07 

lu -0.08 


5300 





130 


47 





1 





72 ±0 


1 1 +0.02 
iJ --0.01 


5300 





240 


19 











29 ±0 


07+0.01 
u '-0.01 


5300 





420 


19 











33 ±0 


A7+0.02 
u ' -0.02 


9500 





130 


35 





2 





67 ±0 


1 1 +0.03 
1 -1-0.04 


9500 





240 


40 











64 ±0 


1 n+0.02 

iu -0.03 


9500 





420 


15 











25 ±0 


06 +0.02 
uu -0.01 


17000 





240 


19 





1 





44 ±0 


1 n+0.03 

±u -0.03 


17000 





420 


8 








n 1 c+o.08 +0.01 
u - iJ -0.05 -0.01 


30000 





420 


5 








n 1 A +0.09 +0.02 
u ' iq: -0.06 -0.02 



Table 2: Values of the reduced cross section. The following quantities are given: 
the values of Q 2 and x at which the cross section is quoted; the number of data 
events, -/V D ata/ the number of expected background events, N BG and the measured 
cross section, with statistical and systematic uncertainties. 
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da/dQ 2 


Q 2 (GeV 2 ) 


da/dQ 2 (pb/ GeV 2 ) 


4tat (%) 


S sys (%) 


5unc (%) 


Si (%) 


s 2 (%) 


s 3 (%) 


280 

OoU 

950 

1700 

3000 

5300 

9500 

17000 

30000 


2.51 ■ 10~ 2 

Z. ( u ■ 1U 

2.13 ■ 10~ 2 
1.22 • 10~ 2 
6.32 ■ 10- 3 
2.88 • lO" 3 
1.24- lO" 3 
2.37- 10" 4 
2.68 • 10~ 5 


±25 

-4-1 8 
±11 

±8.8 
±8.6 
±9.3 
±11 
±18 

+55 
-37 


+21 
-18 
+13 
-16 
+6.0 
-5.9 
+3.4 
-3.6 
+2.0 
-2.6 
+3.5 
-3.7 

+3.3 
-1.1 

+6.7 
-5.5 

+11 
-13 


+20 
-18 
+10 
-13 
+4.0 
-4.7 
+2.5 
-2.9 
+1.5 
-2.5 
+1.2 
-2.1 

+0.9 
-2.8 
+2.6 
-1.2 
+0.8 
-5.4 


+5.1 
-0.6 

-1.3 
-3.4 

+3.4 
-1.7 
+0.0 
-0.3 
+0.7 
+0.0 
-0.7 
+1.3 
-3.0 
+2.9 
-5.3 
+5.9 
-11 
+ 11 


+2.3 
+5.5 

-1.5 
-4.1 

-0.7 
+0.4 

+0.2 
+0.7 
+1.0 
-0.3 
+0.2 
+0.2 

+0.3 
-1.1 

-0.5 
+1.6 
-5.5 
+2.6 


-3.7 
+3.7 
-8.0 
+8.0 
-3.0 
+3.0 
-2.1 
+2.1 

-0.6 
+0.6 
+3.0 
-3.0 
+1.2 
-1.2 
+0.4 
-0.4 
+1.9 
-1.9 


da/dx 


X 


da/dx (pb) 


<Wt (%) 


4ys (%) 


Sunc (%) 


Si (%) 


s 2 (%) 


S 3 (%) 


0.015 
0.032 
0.068 
0.130 
0.240 
0.420 
0.650 


6.08- 10 2 
4.23 • 10 2 
2.76 • 10 2 
1.87- 10 2 
8.08- 10 1 
2.48 • 10 1 
1.28 • 10° 


±18 
±11 
±7.8 
±7.8 
±8.9 
±14 

+137 
-67 


+10 
-9.1 
+5.2 
-7.7 
+1.7 
-2.6 
+2.2 
-1.4 

+2.3 
-2.1 

+5.8 
-4.1 
+15 
-13 


+8.5 
-9.1 
+2.6 
-6.3 
+1.6 
-2.6 
+1.2 
-0.6 
+1.6 
-1.5 
+2.1 
+0.6 
+4.6 
-2.6 


+4.9 
-0.9 
+1.0 
-0.7 
+0.1 
-0.4 
-0.2 
+1.2 

-1.3 
+1.3 
-3.6 
+4.9 
-12 
+13 


+2.5 
-0.3 
+0.9 
+0.1 
+0.1 
+0.0 
-0.8 
+1.0 
-0.2 
+0.4 
+2.2 
-1.8 
+5.8 
-4.1 


-0.1 
+0.1 
-4.3 
+4.3 
-0.3 
+0.3 
+1.0 
-1.0 
-0.8 
+0.8 
-0.2 
+0.2 
+0.8 
-0.8 


da/dy 


y 


da/dy (pb) 


4tat (%) 


Ssys (%) 


5unc (%) 


Si (%) 


s 2 (%) 


5 3 (%) 


0.05 
0.15 
0.27 
0.41 
0.55 
0.69 
0.83 


1.40 • 10 2 
1.18 • 10 2 
8.85- 10 1 
5.63- 10 1 
5.40 • 10 1 
5.07- 10 1 
3.01 ■ 10 1 


±12 
±8.4 
±8.2 
±10 
±11 
±12 
±18 


+8.7 
-1.9 
+4.0 
-3.9 
+2.9 
-4.9 
+3.6 
-2.7 
+3.3 
-3.6 
+3.7 
-5.9 
+12 
-9.9 


+8.6 
-1.8 
+1.8 
-2.2 

+1.5 
-3.3 
+2.8 
-2.6 
+2.7 
-3.0 
+2.5 
-4.7 
+5.8 
-8.3 


+1.1 

-0.4 
+0.6 
+0.3 
-0.8 
-0.3 
-0.1 
+1.0 
-0.1 
-0.3 
-1.9 
+1.6 
-2.3 
+5.7 


-0.4 
+0.9 
+1.0 
+0.5 
-0.2 
-2.4 
+1.9 
+0.3 
+0.3 
-0.7 
-2.3 
+1.4 

-1.7 
+6.7 


+0.1 
-0.1 
-3.3 
+3.3 
-2.5 
+2.5 
-0.8 
+0.8 
+1.8 
-1.8 
+1.8 
-1.8 
+4.6 
-4.6 



Table 3: Values of the differential cross-sections da/dQ 2 , da/dx and da/dy. The 
following quantities are given: the value at which the cross section is quoted; the 
measured cross section; the statistical uncertainty; the total systematic uncertainty; 
the uncorrelated systematic uncertainty and those systematic uncertainties with 
significant (assumed 100%) correlations between cross-section bins. The systematic 
uncertainties considered to be correlated were: the first estimate of the calorimeter 
energy-scale uncertainty (Si); the second such estimate (S 2 ) (see text); and the 
uncertainty in the parton-shower scheme(S?„) . 
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Q 2 (GeV 2 ) 


X 


a 


<Wt (%) 


<*sys (%) 


<Uc (%) 


Si (%) 


82 (%) 


S 3 (%) 


280 





032 





64 


+57 
-39 


+21 
-9.6 


+18 
-9.3 


+5.7 
+0.1 


-1.4 
+7.1 


-1.4 
+1.4 


530 





015 


1 


13 


+43 
-32 


+13 
-23 


+12 
-22 


+5.1 
-1.2 


-3.5 
-6.8 


-2.0 
+2.0 


530 





032 





91 


+39 
-29 


+8.8 
-18 


+3.1 

-15 


-1.1 
-0.2 


+2.4 
-4.5 


-7.9 
+7.9 


530 





068 





83 


+50 
-35 


+17 
-22 


+11 
-17 


-5.7 
-4.5 


-2.5 
+1.6 


-13 
+13 


950 





015 


1 


51 


±26 


+17 
-19 


+16 

-17 


+5.9 
-5.1 


-2.8 
—4.0 


-0.8 
+0.8 


950 





032 





75 


±25 


+11 

-11 


+6.3 
-9.7 


+4.6 
-1.1 


+3.6 
+3.8 


-5.1 
+5.1 


950 





068 


1 


01 


±19 


+12 

-13 


+9.5 
-8.7 


+0.3 
-4.0 


-2.7 
-2.2 


-8.0 
+8.0 


950 





130 





79 


±25 


+14 
— 11 


+8.0 
—4.8 


+4.2 
+3.3 


-1.9 
+2.2 


+9.6 
—9.6 


1700 





032 





87 


±18 


+13 
-5.8 


+12 
-5.6 


+0.7 
+1.9 


+2.6 
+2.0 


-1.4 
+1.4 


1700 





068 





81 


±14 


+2.6 
-5.9 


+1.8 
-5.4 


+0.7 
-1.0 


+0.5 
-1.4 


-1.7 
+1.7 


1700 





130 





75 


±18 


+3.5 
-6.3 


+1.3 
-4.9 


-0.3 

-1.8 


-1.3 
+0.2 


-3.3 
+3.3 


1700 





240 





47 


±24 


+2.7 
-4.8 


+1.1 

-3.0 


-1.2 
-0.3 


-2.9 
+1.3 


-2.1 
+2.1 


3000 





068 





71 


±14 


+6.3 
-1.9 


+4.7 
-0.3 


+2.5 
+1.6 


+2.1 
+1.1 


-1.9 
+1.9 


3000 





130 





58 


±17 


+5.9 
-3.7 


+4.9 
-2.3 


+1.0 
-0.1 


+0.9 
+0.7 


+2.9 
—2.9 


3000 





240 





49 


±19 


+1.7 
-3.0 


+0.9 
-2.4 


-0.5 
-0.6 


+0.5 
-0.8 


-1.4 
+1.4 


3000 





420 
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+49 
-35 


+11 

-9.8 


+8.9 
-9.7 


-0.9 
+5.6 


+2.5 
+0.2 


+1.5 
-1.5 


5300 





068 





53 


±18 


+13 

-15 


+3.6 
-7.5 


-2.1 
+0.1 


-1.7 
+1.1 


+12 
-12 


5300 





130 





72 


±15 


+2.7 
-1.1 


+1.5 
-1.1 


+0.1 
+1.7 


+0.1 
+1.3 


+0.0 
-0.0 


5300 





240 





29 


±23 


+2.8 
-5.1 


+1.7 
-4.9 


+0.2 
+1.6 


+1.5 
-1.2 


-0.5 
+0.5 


5300 





420 





33 


±23 


+6.0 
-6.2 


+2.3 
-2.8 


-1.0 
+2.3 


+2.5 
-3.2 


-4.4 
+4.4 


9500 





130 





67 


±17 


+3.8 
-6.2 


+1.4 
-4.9 


-3.3 
+3.5 


-2.0 
+0.2 


+0.0 
-0.0 


9500 





240 





64 


±16 


+2.4 
-4.7 


+0.4 
-3.6 


-2.3 
+1.1 


+1.1 

-1.1 


+1.8 
-1.8 


9500 





420 





25 


±26 


+6.3 
-5.1 


+2.5 
-0.4 


-3.5 
+3.9 


+4.1 
-3.5 


-0.9 
+0.9 


17000 





240 





44 


±23 


+7.1 
-6.1 


+2.7 
-1.6 


-5.0 
+4.8 


-3.1 
+4.3 


-0.9 
+0.9 


17000 





420 





15 


+50 
-35 


+9.6 
-7.3 


+3.1 
-1.3 


-4.9 
+6.2 


+6.2 
-4.6 


+2.3 
-2.3 


30000 





420 





14 


+68 
-43 


+13 
-13 


+1.8 
-5.4 


-9.4 
+10 


-5.9 
+5.1 


+5.3 
-5.3 



Table 4: Values of the reduced cross section. The following quantities are given: 
the values of Q 2 and x at which the cross section is quoted; the measured cross sec- 
tion; the statistical uncertainty; the total systematic uncertainty; the uncorrelated 
systematic uncertainty and those systematic uncertainties with significant (assumed 
100%) correlations between cross-section bins. The systematic uncertainties con- 
sidered to be correlated were: the first estimate of the calorimeter energy-scale 
uncertainty (Si); the second such estimate (82) (see text); and the uncertainty in 
the parton-shower scheme(S^). 
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Figure 1: Comparison of the final e~p CC data sample (solid points) with the 
expectations of the sum of the signal and ep background Monte Carlo simulations 
(light shaded histogram). The ep background Monte Carlo is shown as the dark 
shaded histogram, (a) the missing transverse momentum, PT,miss, (b) -Pr.miss ex- 
cluding the very forward cells, -Remiss? (c) the variable 5, defined in the text, (d) 
the ratio of missing transverse momentum to total transverse energy, PT.miss/ Et, 
(e) 7/i, (f) the number of good tracks, (g) the Z position of the CTD vertex for the 
high-'jo sample and (h) the Z position of the timing vertex for the Iow-^q sample. 
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Figure 2: The e~p CC DIS cross-sections (a) da/dQ 2 , (b) da/dx and (c) da/dy 
for data (solid points) and the SM expectation evaluated using the ZEUS NLO QCD 
fit, CTEQ5D and MRST(99) PDFs. The insets show the ratios of the measured 
cross sections to the SM expectations evaluated using the ZEUS NLO QCD fit. 
The statistical errors are indicated by the inner error bars (delimited by horizontal 
lines), while the full error bars show the total uncertainty obtained by adding the 
statistical and systematic contributions in quadrature. The shaded band shows the 
uncertainties associated with the PDFs estimated using the ZEUS NLO QCD fit. 
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Figure 3: The reduced cross section, a, as a function of Q 2 , for different fixed 
values of x. The points represent the data, while the expectation of the Standard 
Model evaluated using the ZEUS NLO QCD fit is shown as a line. 
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Figure 4: The reduced cross section, a , as a function of x, for different values 
of Q 2 . The points represent the data, while the expectation of the Standard Model 
evaluated using the ZEUS NLO QCD fit is shown as a solid line. The separate 
contributions of the PDF combinations x(u + c) and (1 — y) 2 x(d + s) are shown by 
the dashed and dotted lines, respectively. 
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Figure 5: The reduced cross section, a, as a function of(l—y) 2 , for different fixed 
values ofx, fore~p (solid points) and e + p (open circles) CC DIS. The expectation 
of the Standard Model evaluated using the ZEUS NLO QCD fit is shown as a solid 
line. The contributions of the PDF combinations x(u + c) and x(u + c) are shown 
by the dashed and dotted lines, respectively. 
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